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Abstract 


This  paper  summarizes  the  status  of  a  direct  construction 
.of  an  asymptotic  representation  of  a  modulating  multiphase 
wavetrain  for  a  class  of  perturbed  KdV  equations.  This  class 
includes  the  KdV-Burgers'  equation.  The  calculations  apply 


on  a  boundary  between  dispersive  and  dissipative  behavior. 

The  construction  proceeds  by  standard  asymptotic  methods.  The 
result  of  the  construction  is  an  invariant  representation  of 
the  reduced  equations  which  permits  their  diagonalization. 

While  mathematically  the  construction  is  incomplete,  care  is 
taken  to  correctly  identify  the  mathematical  status  of  each 
step  in  the  construction.  The  equivalence  of  this  constructive 
approach  with  the  postulated  averaging  of  conservation  laws  is 


established  for  two  phase  waves. 

1 


I .  Introduction 


In  this  paper  we  study  the  initial  value  problem  for  a 
perturbed  Korteweg  de  Vries  (KdV)  equation: 

V  60Dx  +  '2  Uxxx  +  6  f(0'eV  e2uxx"">  -  0 
U(x,t=0)  =  Uin  <  I  »  X  )  . 

Here  e  denotes  the  small  parameter  (0  <  e  <<  1) while 

6  >  0  is  a  constant  of  order  unity.  The  initial  data  U. 

—  ■  xn 

depends  upon  two  spatial  scales ,  a  fast  scale  x/e  and  a  slower 
scale  x  .  The  reason  we  tie  the  fast  scale  in  the  initial  data 
to  the  small  parameter  e  in  the  equation  is  that  the  KdV 
equation  (0  =  0)  in  the  "small  dispersion  limit"  (0  <  e  <  <  1) 
is  known  to  develop  oscillations  whose  spatial  wavelength  is 
0(e).  Actually,  we  will  restrict  our  attention  to  a  special 
class  of  initial  data  -  a  slowly  varying  N-phase  waveform  - 
which  we  will  describe  in  detail  later. 

Our  goal  is  to  construct  a  representation  of  the  solution 
U  =  (x ,t)  which  is  valid  for  small,  but  finite,  e  . 

There  are  two  general  reasons  for  our  interest  in  this 
problem.  The  first  is  technical.  We  desire  to  develop  a  method 
for  the  construction  of  U'  '  which  (i)  is  general  enough 
to  include  external  perturbations  (such  as  dissipation)  of  a 
completely  integrable  system,  (ii)  shows  the  validity  of  our 
prescription  [1,2,3]  of  "averaged  conservation  laws"  for  the 
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description  of  modulating  N  phase  wavetrains,  (iii)  is  a 

sufficiently  standard  mathematical  technique  that  the  construction 

can  provide  a  first  step  toward  a  rigorous  proof  of  the  validity 

of  the  representation.  There  are  cases  (the  sine-Gordon 

equation,  [4,5]  for  example)  where  the  wave  train  is  modu- 

lationally  unstable.  To  study  this  instability  one  must  retain 

additional  terms  in  the  modulation  equations.  The  technique 

we  are  developing  is  sufficiently  standard  that  a  systematic 

retention  of  such  higher  order  terms  can  be  attempted. 

The  second  reason  we  are  interested  in  this  problem  concerns 

the  Korteweg-de  Vries  -  Burgers'  equation  itself.  Consider 

2 

the  special  perturbation  f  *  -  e  Uxx  .  Then  the  equation 
becomes 


2 

0  -  6  UU  +  e  U  -  6  U  =  0 

ut  X  XXX  XX 


2 

with  6=$  e  .  Relax  this  relation  between  5  and  e  for 
the  moment  and  think  of  the  equation  as  depending  upon  two  in¬ 
dependent  small  parameters,  e  and  6  .  When  the  dispersive 
term  is  absent  (e  =0),  the  solution  U  tends,  as  6+0,  to 
a  sequence  of  shocks  whose  speeds  are  fixed  by  the  entropy 
condition.  On  the  other  hand,  when  the  dissipative  term  is 
absent  (6  =  0) ,  the  solution  U  becomes  very  oscillatory  for 
small  £  .  In  this  case  the  weak  limit  of  U'  as  £  +  0  is 
described  by  a  family  of  nonlinear,  hyperbolic  equations  [6, 7, 8, 9]. 
When  both  the  dissipative  and  dispersive  terms  are  present,  but 
small  (0  <  e,  6  <  <  1) ,  shocks  will  describe  the  weak  limit  provided 


SL 


the  dissipation  dominates.  Recently  it  has  been  established 
[  10  ]  that  shocks  will  describe  the  weak  limit  as  e  +  0,  64-0 , 


provided  the  dissipation  parameter  i >  M 


Our  calculations 


(6=8  e  )  describe  the  boundary  between  dissipative  and 
dispersive  behavior.  For  stronger  dissipation,  shocks  apply; 
on  this  boundary,  the  behavior  for  Small  e  is  described  by 
the  equations  derived  in  this  paper;  for  still  weaker  dissipation, 
the  weak  limit  is  described  by  [6, 7, 8, 9]. 

0 

The  method  that  we  are  developing  has  several  main  advantages: 
(i)  It  constructs  a  representation  of  the  wave  which,  since 
it  is  valid  for  finite,  but  small  e,  retains  the  oscillatory 
structure  of  the  wave.  (ii)  It  extends  naturally  to  small 
perturbations  of  integrable  systems. 

This  work  has  two  shortcomings.  The  first  concerns  a 
limitation  of  the  method  itself.  It  demands  a  very  restricted 
class  of  initial  data-slowly  varying  N-phase  waves.  This  means 
that,  for  general  data,  it  is  not  valid  uniformly  in  space; 
rather,  it  applies  only  in  the  vicinity  of  the  "shock  front". 

The  second  shortcoming  concerns  the  work  reported  herein;  not 
the  method  itself.  This  work,  although  systematic,  is  not 
rigorous;  we  view  it  as  the  first  step  toward  a  rigorous  derivation. 
In  the  text  we  will  clearly  , identify  those  points  where  rigor 
is  absent. 

This  work  continues  our  study  of  the  modulation  of  completely 
integrable  wave  trains  [1,2, 3, 4, 5],  The  procedure  of  averaging 
conservation  laws  to  obtain  modulation  equations,  as  well  as  placing 
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these  equations  in  Riemann  invariant  form,  was  initiated  in 
[11,12]  for  the  single  phase  case.  The  constructive  method 
employed  here  follows  single  phase  work  in  [13,14]  and  the 
multi-phase  procedure  of  [15].  In  the  single  phase  case  the 
results  are  not  new,  although  our  arguments  are  somewhat  more 


4 


systematic;  in  the  multiphase  case,  our  results  are  new  because 
of  our  use  of  the  completely  integrable  exact  theory.  In  the 
single  phase  case,  dissipative  perturbations  of  nonlinear. 


dispersive  wave  trains  have  been  analysed  previously  [12,16,17]. 


The  (weak)  zero  dispersion  limit  of  the  KdV  equation  is  studied 
in  [8,9,10]. 
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The  last  section  of  this  paper  is  motivated  by  the 
rigorous  work  of  [19,20,10]. 

Finally,  I  acknowledge  many  conversations  with  R.  DiPerna, 


H.  Flaschka,  M.  G.  Forest,  C.  D.  Levermore,  H.  McKean,  G.  Papanicolao 
and  S.  Venekidas.  My  work  has  certainly  benefited  from  each  of 
these  Interactions. 


II.  Definition  of  the  Problem 


Consider  the  initial  value  problem 

°t  -  6  DDx  +  **  °xxx  +  6  f(U'  «V  ^uxx'-">  - 

.  .  *+■ 

U  (x,t  =  0)  =  wPJ>  (  |  ,•  X  (x)  )  . 


0  , 

(II.la,b) 


Here  0  <  e  <  <  1  and  $  =  0(1)  .  The  initial  data  is  a 
slowly  modulating  N-phase  waveform  for  the  KdV  equation  which 
we  now  describe. 


The  KdV  equation  has  a  family  of  exact  solutions  of  the 

form 


»©*  (x,t)  »©Ux,t) 

U  (x,  t)  —  Wjj(  - ^  t  •  •  •  >  ^  #  X^, . . . ,  X^  )  II. 2a) 

where  the  N  "phases"  ^,...,>6^  depend  linearly  upon 
x  and  t  ,  • 


»Oij  (x,t)  5  !Cj  x  “  «j  t  .  (II. 2b) 

The  waveform  WN  is  2  it  periodic  in  each  argument 

0.  =  KH/e  .  The  family  of  "N  phase  waves"  (II. 2a)  is 
J  D  4 

indexed  by  2N+1  parameters  X  =  (XQ,  X1#...fX2N)  ,  which 
fix  the  N  spatial  wavenumbers  ic  =  (k^,  . . .  ,kn)  and  the  N 
temporal  frequencies  u  =  (<d^,...,un)  .  The  explicit  formulas 
for  k  *  k (  X  )  and  u  ■  w(t),  as  well  as  more  details  about 
these  N  phase  waves,  will  be  given  later.  For  now,  we  only 
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remark  that  in  the  single  phase  (N=l)  case,  (II. 2)  is  a  2tt 

periodic  traveling  wave  solution  of  the  KdV  equatiom. 

The  initial  data  (II. lb)  is  described  in  terms  of  an 

N  phase  wave  W  as  follows:  Prescribe  smooth  functions 
-*■  ■ 

X  =  X (x)  in  terms  of  which 


where 


(II. 3a,b) 


Vj 


Kj  (X  (x)  ) 


X (x)  prescribed 


Our  problem  is  to  construct  a  representation  of 
U  -  Ue  which  is  valid  for  small  e  .  We  use  a  standard 


asymptotic  method  which  begins  with  an  ansatz: 
fcf(x,t,  *©£(x,t)  1 

ovu  Ft- . ^-r-»  H  • 


(II. 4) 


Here  the  (x,t)  dependence  of  *©^  *  («©^,...^©^)  is  to  be 


determined.  We  define 


•ce(x,t)  =  3x»0.e(x,t) 


we(x,t)  s  3£©£(x,t)  , 


and  note  that 


3t  °e  ^  +  3t)uE 


a  uc  (JL.  1  +  a  )ue  , 
x  e  x 


(II. 5) 


where  V  U  =(^11  ,...,  -g^U  )  “  e  ^ae  u  •  •  •  •  •  3©  u  ^*  of 

in  IN 

course,  definition  (II. 5)  implies  N  consistency  (or  integrability) 

conditions  for  6  . : 

3 


H' 


3.  Ke  =  3  Ue  . 
t  X 


In  terms  of  this  ansatz,  equation  (Il.la)  becomes 


ue'Vue  -  6Ue  Ke.VUe  +  (ice*V)3ueJ 


(II. 6) 


[°t 


.  n>  ^  1 


4 


+  lur  +  6U"U^  +  3l(Ke-«V)‘  U^+dc^.V)  (ic^*V)Ufc|+3f  (Ue,  (Ke*v+e3 


x)lle,...)] 


[■ 


+  e | Ke*VUe  +  3  kE  *vue  +  Kc  «VUe|  +  e2 

1  XX  XX  XX  1 


:"1  +  e2  [ue  1  = 
J  xxx 


=  0  . 


In  this  equation,  expansions  of  the  form 


UG  ^  W  +  e  U<1J  +  e2  U*2>  +  ... 


k  *v  k  +  e  IC 


(1>  +  s2  k<2>  +  ... 


(XX. 7) 


(0  *v*  w  +  e 


»(1>  +  e2  u<2>  +  ... 


lead  to  the  following  sequence  of  problems: 


(He”1) :  (wV)  W  -  6  W(k*V)W  +  (tc*V)3  w  =  0 

O(e^)  :  L  U(:5)  +  P(^  »  0  ,  j  -  0,1,2, 


(II. 8 
(II. 8 


■  v*  O  -  * 

Li'A', 


9 


45555555 55555 


,.«v. %.* v *v  v\. *: .  />* 


•  • 


where  the  linear  operator  L  is  given  by 


L  =  oj* V  -6  (k  *V)  W  +  (k  *V)  . 


The  inhomoge neity  Fl  '  is  given  by 


(II. 9) 


=  [wfc  -  6WWx+3  (U*V)2Wx  +  (k  •  V)  (  <x*  V) wj  +6  f  (W,  (K*V)W,...)j 


+  [u)(1)*VW  -  6  W  <{1)*VW  +  3  (k  •  V) 2  (K(1)*V)  W]  ,  (II. 10 


with  similar,  but  more  complicated,  formulas  for  F1-5'  ,  j  >_  1. 

The  sequence  of  problems  (II. 8)  must  now  be  studied;  however, 
we  first  summarize  some  background  material  from  the  theory 
of  the  inverse  spectral  representation  of  KdV.  This  material 
will  be  used  throughout  the  paper. 


ZJ*. 
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III.  Inverse  Spectral  Theory  for  the  KdV  Equation 


The  background  material  in  this  section  may  be  found  in 
[1,21]  which  contain  references  to  the  original  literature. 
Let  q  =  q(y,T)  satisfy  the  KdV  equation, 


q  =  6  qq  -  q 

M'iy  ^yyy 


(III.l) 


By  considering  the  "Lax  representation"  of  this  equation, 
one  realizes  that  it  arises  as  the  integrability  condition 
for  the  linear  system 


[-3  +  2(q3  +  3  q)]  Y  =  4X  3  ¥ 

yyy  *  y  y^1  y 


3X  f  «  -  2qyY  +  2(q+2X)3yY  . 


(III. 2a, b) 


The  function  ¥  is  known  as  a  " squared  eigenfunction" . 
Equation  (III. 2a)  is  an  eigenvalue  problem  for  the  squared 
eigenfunction  ¥ ;  equation  (III. 2b)  defines  its  time  flow. 

The  pair  (III. 2a, b)  is  compatible  since  q  satisfies  KdV. 

The  pair  of  equations  (III. 2)  is  fundamental  in  the  theory 
of  the  KdV  equation.  Here  we  use  the  pair  (i)  to  generate 
an  infinite  family  of  conservation  laws,  (ii)  to  provide  a 
representation  of  the  N-phase  wave  solutions  of  (II. 8a) , (iii) 
to  provide  formulas  for  the  fluxes  and  densities  of  the 
conservation  laws  in  terms  of  the  N-phase  waves,  and  (iv)  to 
compute  necessary  averages. 
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III. A.  Fundamental  Conservation  Law 


Let  q  solve  KdV  and 


solve  the  pair  (III. 2).  Then 


it  immediately  follows  that  Y  satisfies  the  conservation 


3Tm  +  ay  [6(q-2X)H'  -  23^]  =  0 


(III  .3) 


This  fundamental  conservation  law  generates  an  infinite  family 
as  follows:  One  seeks  a  solution  of  (III. 2a)  which  has 
the  asymptotic  behavior 


Y(y,T;A)  =  +  X]  Y.  (y,T)  (2A)"j“1/2  as  X  -  «.  (III. 4) 


This  ansatz  in  (III. 2a)  leads  to  a  recursion  relation  for  the 
coefficients  'Pj(y,T), 


8y  Vl  =  7  8yyy  +  *>8y  +  8y  1 


(III. 5) 


%  -  1  • 


Recursion  relation  (III. 5) ,  together  with  certain  choices  of 
integration  constants,  provides  formulas  for  in  terms  of 

q  and  its  derivatives.  We  list  the  first  few: 


7  <3<j2-v 

5  3  .  5  n2 

J  q  *  +  T  qy 


^y^7  ^y  T  ^yyy^ 


35  4  35  2  7  2  L  _ 

“5"  *1  +  ~r~  qq  +  s-  q  +3 

8  4  v  8  nw 


+  I  <M, 


-  7-  q  q 
4  y  yy 


»,  [-  * 


c  2 

5q  qx 


By  inserting  the  asymptotic  behavior  (III. 4)  into  the  fundamental 
conservation  laws  (III. 3),  we  obtain  a  family  of  conservation 
laws  for  the  KdV  equation: 


W  +  ySWj  -  -  “yy^j1  =  0  -  i  *  1.2 


(III. 7) 


III.  B.  Solutions  of  the  Adjoint  Linearized  KdV  Equation 

The  squared  eigenfunction  Y  generates  solutions  of  the 
adjoint  linearized  KdV  equation.  Let  q  and  Q  =  q+6q  denote 
two  solutions  of  the  KdV  equation.  Then,  as  their  difference 


6q  goes  to  zero,  it  satisfies  the  linear  equation 


3T  «q  -  6  3y  q6q  +  3yyy  5q  =  0 
with  formal  adjoint 

#t  ’  -  s^sy  '  *  3yyy  'K  =  0 


(III. 8) 


(III. 9) 


Let  q  denote  a  solution  of  KdV,  and  Y  denote  a  solution 
of  the  pair  (III. 2).  Then  it  immediately  follows  that  Y 
satisfies  the  adjoint  equation  (III. 9) ;  thus,  the  coefficients 
Yj  are  also  solutions  of  the  adjoint  equation. 
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III.  C.  The  v  Representation  of  KdV  Waves 

The  squared  eigenfunction  system  (III. 2a, b)  can  be  used 
to  generate  representations  of  N-phase  watetrains.  This 
construction  begins  with  the  observation  that  system  (III. 2a, b) 
admits  a  first  integral  [22] : 


1  YY  -  1  V2  -  (q-X)  Y2  =  R2(X) 

2  yy  4  y 


(III. 10) 


To  construct  N-phase  waves,  we  fix  (2N+1)  real  constants 
X  =  (X^X^...,!^)  , 


X0  <  X1  <  ***  <  X2N  ' 


(III. 11) 


and  demand  that  the  first  integral  R  (X)  be  the  polynomial 


R  lX)  =  n  (X-X  )  . 
k=0  k 

This  situation  is  achieved  by  seeking  a  solution  Y 


(III. 12) 


of  the 


squared  eigenfunction  system  which  is  polynomial  in  X  of 


degree  N  : 


(y,T;X)  =  n  fx-v.(y,T)l  . 

j=i  l  3  J 


(III. 13) 


Inserting  this  ansatz  into  (III. 2)  leads  to  the  "y-representation 
of  KdV  waves  with  N  degrees  of  freedom"  : 


••  r* 


Formulas  (III.14)  summarize  the  "y-representation  of 
KdV  waves  with  N-degrees  of  freedom".  We  now  show  that  these 
exact  solutions  of  KdV  are  "N  phase  wavetrains"?  that  is, 
they  are  solutions  which  are  quasi-per iodic  in  space  and 
time  which  depend  upon  N  phases. 


III.  D.  The  Representation  of  KdV  Waves 

The  wave  qN  admits  an  equivalent  representation  which 
results  once  an  Abel  transformation  is  used  to  integrate 
the  y  equations  (III.14d,e).  On  the  Riemann  surface 

/m - 

R  “  [A,R(A)  =  /  n  ( X— X. ) 3  ,  one  fixes  a  canonical  set 

0  K 

of  a-b  cycles  (Figure  1).  On  this  Riemann  surface,  we  introduce 
the  following  objects:  (i)  a  basis  of  holomorphic  differentials, 


Figure  1 


♦i 5 


o-l  dX 

cij  x  rTTT  * 


(III. 15a) 


normalized  by  the  condition 


f.  *j  *  4ij  ’ 


(III. 15b) 


(ii)  A  symmetric  matrix  with  positive  definite  imaginary  part. 


3  H  L  *3  ' 


(III.lSc) 


(iii)  Two  differentials  of  the  second  kind 


dX 

RTXT  ' 


°15  [-J  +  §°3  xJ*1] 


r  n  -i 

,  .N+l  .1*1*.  \  1  K  ! j-1  dX 

n2s  - 1  *  i A  x  ^i*  3 


(III. 16a, b) 


where  the  coefficient  A  = 


X^  ,  and  the  coefficients 


and  Ej  are  uniquely  determined  by  the  normalization  conditions 

o  fl.  *  0 
Ibi  1 

These  quantities  form  the  ingredients  for  a  change  of 
variables  from  (p^,...,y^)  -*•  (0^,...,0jj)  t 
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e, 


(b”1) 


ij 


N  1 

2  f 

k=l 


(III. 17  ) 


If  lij(y/T)  satisfies  the  differential  equations  (III.14d,e)  , 
then  the  new  variables  0j(y,T)  satisfy 


9y  ‘  "5 


3t  ej  "  Wj  ' 


(II. 18a, b) 


where  the  constants  k.  and  oi.  are  defined  in  terms  of 

4  J  D 

X  : 


-  kj  <*>  -  - 1  °i 


0>J  *  U)j  (X)  »  -  12 


L  °2  • 


(II. 19a, b) 


Using  theta  functions,  one  can  invert  transformation  (III. 17) 
and  give  formulas  for  (and  therefore  q^)  in  terms 

of  (0  j,  •  •  • » 0jj )  s 


qN  -  A  +  T  -  23^  log  |*e»(z  (0)  ;B)^  , 

M 

where  the  theta  function  is  defined  (for  z  e  C  )  by 


(III. 20a) 


•6>(z;B)  =  ^  i  exp{tti[2  (m,z)  +  (m,Bm) )}  , 

me?* 


(III. 20b) 


(III.20c) 


and  z (6)  denotes  the  linear  map 


zj(0)  -  I  el 


where  is  the  j—  column  of  the  period  matrix  B 


Finally,  the  constant  T  is  given  by 


N 

r  E  -2  S  L x 


(III. 20d) 


and  d  denotes  a  real  constant  which  plays  no  role  in  the 
following. 

In  this  manner  the  theory  of  the  exact  KdV  equation  has 
generated  a  solution  (III. 20a)  in  the  form 


qN  "  qN(0l'***'0N;X) 


where  the  (y,x)  dependence  enters  only  through  the  phases. 


©j(y#t)  *  K^y  +  +  6^  , 

and  where  the  wave  form  is  2n  periodic  in  each  individual  phase. 
Thus,  the  constants  and  are  enterpreted  physically 

as  spatial  wave  numbers  and  temporal  frequencies. 


This  completes  the  summary  of  that  material  from  the  in¬ 
verse  spectral  representation  of  exact  KdV  waves  which  we  need 
for  this  paper.  Now  we  return  to  the  analysis  of  the  sequence  of 
problems  (II. 8). 


ks-wwaiwaasi 


IV.  The  Leading  Order  (0(e  ))  Problem. 


In  this  section  we  construct  solutions  of  (II. 8a).  Fix 

2N+1  real  constants  X  =  (X^  <  X^  <  . . .  <  X2N)  ,  and 

construct  the  B  matrix  and  the  differentials  .  as  in 

J  -* 

(111,15)  and  (III. 16),  Define  the  wave  vector  k-k(X)  and 
the  freqeuncy  vector  oj  =  w(X)  by 


:i(X)  =  -  | 


<o^(X)  *  “12  o  ^2 


(IV. la,b) 


For  this  k=ic(X)  and  w=&)(X)  ,  we  seek  W:T  (N-Torus)  -*■  ]R 
which  satisfies  the  0(e-1)  problem: 


(U'V)H  -  6W(ic*V)W  +  (<•  V)  W  =  0 


(IV. 2) 


Using  material  in  Section  III,  we  find  solutions  of  (IV. 2) : 


W  =  W(0;O,X)  =  qN(0-0;X) 


(IV. 3) 


Equivalently,  we  may  use  the  y-representation: 


W  =  W(0;0,X)  ■  A  -  2 


VjO-OjX)  , 


(IV. 4) 


where  the  y  variables  satisfy 


R(y4> 


(k*  V)  y  .  =  -2i  -jj - J - 

3  Jj  <>vV 


(IV. 5a) 


(u-y)u,  =  ->‘2i  2A  -  2 


SR(y.) 

n  (Vi,  -  VI.) 

<  V-i  J  1 


(IV.  5b) 


In  this  manner  we  generate  a  3N+1  real  parameter  family  of 

M 

solutions  W:T  -*■  3R  of  (IV. 2).  N  of  the  parameters, 

o  o  o 

0= (0^. . . ,©N)  ,  are  trivial  in  that  they  merely  center  the  N 

phases  0= (0^. . . ,©N)  .  The  remaining  2N+1  parameters 

A  =  carry  qualitative  information  about  the 

wave.  For  example,  they  determine  the  wave  numbers  k-  (ic^  , . . .  ,*N) 
and  the  frequencies  (w^ , . . .  ,n>N)  by  (IV.l).  In  addition, 
they  determine  the  mean  of  the  wave  W: 


<VJ>  = 


Sr  f  w 

!ir)N  J  N 


(0)  d  0  . 


(IV. 6) 


Remark (i)  Physically,  it  seems  more  natural  to  co-ordinatize 

the  N  phase  waves  by  the  N  spatial  wave  numbers  k,  N  temporal 

frequencies  w,  and  mean  <W>  rather  than  by  the  2N+1 

*► 

parameters  A  .  However,  the  X  co-ordinates  are  better  understood 
mathematically , 


Remark  (ii)  For  the  single  phase  (N=l)  case,  equation 
(IV. 2)  is  an  ordinary  differential  equation  of  third  order  which 
is  easy  to  analyse.  One  quickly  shows  that  all  of  its  2n  periodic 
solutions  belong  to  the  4  parameter  family  (IV. 3)  with  parameters 


2 


I 

*  • 

j  «.  » 

! 

*  Remark  (iii)  For  N  >  1,  I  suspect  that  all  solutions  of 

(IV. 2)  which  genuinely  depend  upon  all  N  phases  (|^  t  for  any  j) 

belong  to  family  (IV. 3).  Assume  there  exists  W:TN  ]R  which 

' 

|  solves  (IV. 2),  depends  genuinely  upon  all  N  phases,  and  does  not 

I  belong  to  family  (IV.  3).  Then  W(y,x)  =  W (tc^y+w^x , . . . , K^y+co^x )  is  a 

solution  of  the  KdV  equation  which  is  quasi-periodic  in  y  and  x 
|  with  exactly  N  frequencies,  and  yet  is  not  an  "N-gap  potential" 

for  inverse  spectral  theory.  I  think  that  no  such  solution  exists. 

i  .  Remark  (iv)  Notice  that  in  this  theory  it  is  easy  to 

shut  off  a  phase ,  say  0 .  .  One  seeks  a  W  which  is  independent 

A 

of  0..  If  k.  and  a),  go  to  infinity,  this  situation  is 
3  3  3 

forced  upon  us.  On  the  other  hand,  if  and  ok  vanish 

(a  soliton  limit) ,  3/30^  is  removed  from  equation  (IV. 2) . 

In  any  case,  such  situations  must  be  understood  before  a 
modulation  theory  sufficiently  general  to  create  and  destroy 
phases  can  be  developed. 
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V.  Solvability  Theory  for  the  Linear  Problems 

Fix  W,  a  solution  of  (IV. 2)  in  family  (IV. 3).  All 
of  the  (He-5),  j  0,  problems  are  of  the  form 

LU  +  F  =  0  ,  (V.l) 


with  a  prescribed  inhomogeneity  F:  TN  ->  ]R  .  Here  the  linear 
operator  L  is  defined  in  terms  of  W  by 

LU  =  (wV)U  -  6  ( K  •  V )  W  U+  (k  •  V)  3  U  (V.2) 

(We  work  in  the  Hilbert  space  of  functions  U:TN  -*■  IR  which 

N 

are  square  integrable  over  the  torus  T  . ) 

For  the  solvability  theory  of  (V.l) ,  we  need  to  under¬ 
stand  n  (L)  ,  the  null  space  of  L,  and  n  (L^~)  ,  the  null  space 
of  the  adjoint  of  L.  Here  the  formal  adjoint  is  given  by 

L^V  -  -  (oj*  V)  V  +  6W  (k  •  V)  V  -(<*V)3V  .  (V.  3) 

We  have  the  following  fact  concerning  the  null  space  n (L^) : 


Theorem  VI:  (a)  ¥(N)  =  n  (X-y.)  e  n(L+)  V  X 

j=l  3 


(b) 


^(N)  1  _-i 

TSror  i(2x) 


e  T1(L^)V  X 


(V.4a,b,c) 


(c)  e  n(L')  V  j  =  0,1,... 


23 


The  proof  of  this  theorem  follows  immediately  from  the  material 
around  (III. 14),  together  with  the  fact  that  Uj(y,t)  depends 
upon  (y,t)  only  through  the  phases  Gj(y,t)  =  K^y+u^t. 


Theorem  V.2;  Formulas  (V.4)  generate  only  (N+l)  linearly 
independent  members  of  n  (L-^)  .  These  may  be  represented  as 


N  w 

(a)  {an,a1#...,a  }  ,  where  n  *  /,  a. A3  , 

u  N  3=1  3  3^o  J 


(b)  (v0,ylf...fyN} 


Theorem  V.3:  (a) 


(V. 5a,b) 


e  n(L)  V  j=l , 2 , . . .  ,N 


(V.  fc&,b) 


(b)  1-6  (K.VW  )  W  e  n (L) 


Using  formulas  (V.5),  we  have  (N+l)  solvability  conditions 
which  are  necessary  if  (V.l)  is  to  have  a  solution.  For 
N=1 ,  simple  analysis  of  the  ordinary  differential  equation  shows 
that  these  are  actually  necessary  and  sufficient.  We  have 


Theorem  V.4  Let  N  =  1  .  Then 


(a)  n(L^*)  =  span  (1,W) 


(b)  n(L)  =  span  {W_,  1-6  k  W  } 

6  to 


(V.  7) 


^wVlVl’.  A  .'«  -■“  /* ... 


,>A-w  v  v.v  v,:.  - 


Thus,  in  the  N=1  case,  both  p (L)  and  p (L^)  have 


dimension  2 . 


Remark  (i)  For  N  >  1,  we  suspect  that  p (L)  and 

n(Lt)  have  dimension  N+l.  If  so,  p(L"^")  =  span  {¥g, . . . ,¥^1  . 

We  have  not  succeeded  in  proving  this.  For  now,  the  solvability 
conditions 


(VyF)  =  0  ,  for  j  =  0,1, ...  ,N 


(V.  8 ) 


I 

t 


are  necessary  for  N  >  1;  necessary  and  sufficient  for  N  =  1. 


VI  The  0(e)  Problem 

Armed  with  this  solvability  theory,  we  return  to  the 
sequence  of  linear  problems  (II. 8b).  Explicitly,  the  O(e^) 


problem  is 


L  U<°>  ♦  *<°>  +  -  0  . 


J'<0)  -  [VL  -  6WW  +  3  |(k-7)2W.  +  (k-V)(k  •' 

“  *  I  *  X 


+  ef(w,K‘Vw,...)]  , 


■  [(u)(1^*V)W  -  6W  k^*VW  +  3(k*V)2(k^  *V)W]  . 


In  the  source,  one  does  not  have  to  worry  about  F 
because  for  this  part  of  the  source,  we  have  an  explicit 
solution.  Recall  that  W  satisfies 


3(0) 


(VI, la) 


(VI, lb) 


(rV)W  -  6W(k •  V)W  +  (k * V) ^ W  «  0. 


Define  2N  functions  on  the  torus  Tw  by  . 

y'V  =  IS 

x  *  8V 

(K-j)  gu 

X  3  =  fjjr  ,  i  •  1,2,..  N. 

Then,  by  differentiating  the  W  equation,  one  finds 

K)  a 

LX  3  +  ^  W  "  0 

V 

LX  J  -  6W  ^  W  +  3(k-V)2  ^  WO 

J  J 


(VI, 2) 


(VI, 3) 


Using  (VI, 3)  we  may  write  U 


WV»V  V  .  »  «•  »  •  «  »  •  »•  .  •  .  •  .  •  <  *•*.*, 


as  follows: 


:  *  _  L  ^  _  *  •  »  _  ^  -  ”«'«w  ••«••• 

tA*.  »%•«%*  ■#**---  V3-%A%*  '••  •  •  •  V»%  »  •  «  v  *«. 


5 

1 

2 


•  • 
•  * 


where 


l  [.<*>  x'“J>  4  «j1)  x“  J>3  4  &(0> 

J-l 


L  V(0>  +  ?<°>  .  0 


'v(O) 

Thus,  the  solvability  theory  need  only  treat  p  . 


(VI, 4a) 


(VI, 4b) 


I 


Theorem  (V,2)  provides  N  +  l  necessary  conditions  for 


solvability. 


(Tj  ,  £(0))  ■  0  for  J  -  (0,1,. ..,N). 


(VI, 5) 


In  this  manner,  we  arrive  at  the  modulation  equations  which 
must  be  satisfied: 


Vj  ■  Vj  ’  J-l. 2.  ••••>< 


(VI, 6) 


Vj.  Wt  -  6WWX  +  3  [<k-V)2Wx  +  (K-V)(Kjt.V)W]  +  8f(W,K.VB,..)j 

*  0,J-0,1,...,N 


Equations  (VI, 6)  are  the  main  result  of  this  section. 

They  provide  a  system  of  2N+1  first  order  partial  differential 
equations  for  the  2N+1  parameters  t  *=  (XQ, . . .  »*2N)  •  The 
first  N  equations  result  from  consistency  of  the  ansatz; 
the  last  N+l  from  necessary  solvability  conditions. 

Equations  (VI, 6)  provide  a  closed  system  which  depends  only 
upon  t,  t  ,  In  the  next  section,  we  place  this  system 

in  manageable  form. 

If  the  null  space  is  exactly  N+l  dimensional  (as  we 
know  for  N=1  and  suspect  for  N  >  1),  these  solvability 


■torsi 
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conditions  are  necessary  and  sufficient  to  ensure  the 
existence  of  :  TN  ]R  .  At  this  stage,  one  can 

proceed  in  two  different  directions,  (i)  One  can  continue 
to  generate  higher  order  terms  in  the  expansion  of 

Ue  The  corrections  to  the  frequencies  ajjC  ,  the  wave 
numbers  KjE  ,  and  the  mean  <Ue>  will  provide  sufficient 
freedom  to  ensure  solvability  at  each  stage.  In  this  manner, 
a  formal  asymptotic  expansion  of  the  form 
Ue  +  e^U^^  +  ...  can  be  constructed. 

(ii)  Alternatively,  one  can  truncate  the  expansion  at 
WN  ;  t(x,t))  where  the  t(x,t)  satisfy  the  modulation 

equations  (VI, 6)  and  prove  that  U e  -  WN  is  0(e)  with  some 
uniformity.  The  second  direction  is  the  most  important. 


28 


VII.  Connection  between  the  Modulation  Equation  and 


Averaged  Conservation  Laws 


The  modulation  equations 


*tKJ  "  Sx“j  ' 


]  a  If 2/ ... (N 


(VII. la, b) 


(H* . ,  W+-6WW  +3t(<-V)2W  +(<*7)  (k  -V)W]+Bf  (W,<-Vv.',  ...))  =  0 

J  W  X  X  X 

3  s  0,1,2,... ,N, 

although  a  closed  system  for  £(x,t),  appear  to  be  a 
complicated  system  of  nonlinear  partial  differential  equations. 
In  13]  we  show  that,  even  in  the  presence  of  an  external 
perturbation  f,  these  modulation  equations  are  actually  very 
tractable  provided  (VII. lb)  can  be  replaced  by  N+l  averaged 
conservation  laws.  In  this  section  we  derive  the  validity  of 
this  replacement  for  N  =  1,2.  (This  is  sufficient  to  treat 
2  phase  waves.) 


VII. A.  Averaged  Conservation  Laws 


One  approach  to  deriving  modulation  equations  is  to 
average  conservation  laws.  One  assumes  that  the  exact  equation. 


Ut  -  6UUX  +  xxx  +  8f  (U,cUx,c2Uxx,...)  -  0, 


(VII. 1) 


has  a  solution  of  the  form 


eN<x,t) 


;  t(x,t)j  +  0(e)  ,  (VII. 2) 


%"  +..ML,.+,m  ». 


-.  *l"  *.  *:  v  v« 


which  is  2 it  periodic  in  each  0^. 
conservation  laws  of  the  form 


In  addition,  it  has 


3tT(U,eUx,...)  +  3xX(U,eux, ...  )  +  BG(U,eUx  ,...)  =  0  ,  (VII. 3 


for  any  solution  of  (VII. 1).  In  particular,  evaluating  on 
solution  (VII. 2),  conservation  lav;  (VII. 3)  takes  the 
form  (3t  ■+  i-  w*v  +  3t  ,  etc.) . 


|  [(<*>. V)  T  +  (k*V)X  3 


(VII. 4 


+  T(WN,  (k-V)Wn,...)  +  X(Wn,(k.V)Wn,...)+BG(VJn>«.V)Wn,.. 


+  . . .  *  0 

N  —1 

When  averaged  over  T  ,  the  0(e  )  term  averages  to  zero,  and  one 

is  left  with  the  averaged  conservation  law 


VT(WN'  (k*V)Wn,...)>  +  3x<X(Wn,  (k.V)Wn,...)> 


+  B<G(Wn,  (<»V)Wn  ,...)>  «=  0 


(VII, 


This  is  one  equation  among  2N+1  unknowns. 

Each  KdV  conservation  law  will  lead  to  an  averaged 
conservation  law  of  the  form  (VII. 5).  Indeed,  consider 
j  KdV  density  as  generated  by  the  recursion  relation  (III. 5), 

ki 

=  Yj(U, e3xU,  ...,  (e3x)  ^U)  .  Here  kj  is  the  order  of 

the  highest  derivative  of  U  in  Yj  .  We  compute: 


30 


Vj  '  Jo  TJ>1  9t(eV‘u 

■  Jo  ’M'£Vl'^-t2V  •••» 

k . 

■  -Vj  - s  £  £<u'ev--->  > 


where  7.  #  denotes  a  partial  derivative.  Thus, 
j  •  * 

we  obtain  the  perturbed  conservation  law 

k . 

Mi  +  Mi  =  “B  l  V.  .(eS)8,  f(U,MJ  .  ...)  , 

*-  j  *  J  1=0  J  *  *  x  x 

with  the  right  hand  side  explicitly  given  in  terms  of  the  j 

density  and  the  external  perturbation  f.  Evaluating  this 

conservation  law  on  the  wave  form  (VII. 2),  and  averaging 

N 

over  the  torus  T  yields  the  perturbed  avaraged  conserva¬ 
tion  law 


th 


where 


M'V  +  V6(wrj  “  'rj+i>>  +  ’ 

kj 

<G.>  =  l  (f.  9,(*‘7)A  f(W,K-VW,  ...)), 
3  i=0  j 


(VII. 6a) 


(VII. 6b) 


and  the  densities  are  evaluated  on  the  N  phase  wave  form  W. 
For  later  use,  we  list  the  first  three: 


3t<W>  -  3X<3VT>  +  3<f  (W,  (k«V)W,  . .  .)>  -  0 


M 


,2  .  3 


(VII.a,b,c 


3t<^->  -  3X<2W  +I  (k-7W)‘>  +  3<W  f(W ,k«VW,...)>  =  0 

at<w3+  |  (k«vw)2>  -  ax<|  w4  + 12w(x-vw)2+  |((k«v)2w)2 

+  8<3W2f  +  (tc  •  V) W  (k*7)f>  -  0 


jl 


VII. B.  Connection  between  Averaged  Conservation  Laws 
and  Null  Space  Modulation  Equations 


Consider  the  null  space  equations  (VI. 5), 

(yjfP(0))  *  0  ,  J  =  0,1,2, ...,N, 

(VII. 8a, b) 

P(0)  =  Wt-  6WWx+  3((k-7)2Wx+  (K.V)  (kx-V)W)  +  Bf(W,  (K«V)W,...) 

We  show  the  first  three  of  these  are  averaged  conservation  laws; 
indeed,  they  become  (VII. 6)  for  j  *  1,2,3.  Since  Yq  =  1,  the 
first  is  immediate: 

0  =  (Y0,F(0))  =  (1,  F(0)) 

a  3.  <W>  -  3  <3W2>  +  8<f> 

t  X 

The  next  is  almost  as  simple: 

0  =  ('P1,F(0))  =  (W,F(0)) 

■  (W,W<_)  -  (W,6WW  )  +  3  f  Vf,  (k-V)2W  +  (k»?)  (  k  •  V)  W)  +  8  ( W,  f ) 

*  1  at<w2>  “  2  ax<W3>  -  3(k»Vw,3x(x«V)w)  +  8  (W, f ) 

*  at<T‘>  “  ax<2W3  +  |  (ic*VW)2>  +8<W  f>. 

The  verification  that  the  third  null  space  equation  yields 
the  third  averaged  conservation  law  is  more  tedious. 

It  uses  an  extra  ingredient  —  atKj  55 
tion  begins  as  follows: 


3  w.  .  This  calcula- 

X  J 


0  -  2  (^2#  F  )  -  (3w  -  (  F'w#> 

-  (3W2  -  (  k  •  V )  2  W»  Wt-3(W2)X+  3(k*V)2V^+ 3(ic*V)  (kx«V)  W  +  Bf) 

»  (3W2#Wt)  -  (U*V)2W,Wt)+  9(W2,  (W2)x)  +  3((k-V)2W,  (W2)x) 

+  9(W2,  (k*V)2Wx)  -  3((k*V)2W,  (ic*7)2Wx) 

+  9(W2,  (<*V)  (kx-7)W)  -  3((k*V)2W,  (k*7)  (kx*V)W) 

+  33 (W2,f)  -3 ( (<• V) 2W/ f ) 

=  ^t<W3>  +  (  (k*7)W,  (k  *V)WtJ  -  |  3x<W4>  -  3((k-V)W,  (ic*V)  (W2)x) 

-  9((k*V)W2,  (k*V)wx)  -  3<(k-V)2W,  (k*V)2Wx)  -  9  ( (ic  *7)  W2 ,  V)  W) 

-  3<(k-7)2W,  ik*7)  (k  *7)W)  +  B<  f  [  3  W2  -  (k*7)2W]  > 

2 

*  3t<w3>  +  >  -  ((K.7)w,  (<ct*V)w)  -  3X<  I  w4> 

-  3((tc«7)2V,  (k*7)2wJ  -  *3((ic*7)2W,  (k-7)  (k«7)w) 

*■»  X 

-  3((K*V)Wr(K*V)  (W2)x)  -  9({K*V)WX,  (K’V)W2) 

-  9((k*7)w2,  (kx*V)w)  +  3<f  [3W2-(k*V)2w]  > 


Continuing, 


o  «  3  <W  +  i  (ic»VW)  >  -  3  <i  vr>  +  P<f  [3W^-(k.V)^W]> 

t  &  X  4 

-  (<K*V)W,  Ut*V)w)  -  3X<|  [  ( ic  •  V )  2jWl  2> 

+  3((<*7)2W,  (k*V)  (kx«7)w) 

-  3  <12W(k*VW)2>  +  12  (w,  (ic*VW)  2)  +  24  (w,  t  ( ic  -  7)  W]  ic-VW  ) 

AX  X 

+  24(W,  [  (  *c  *  V )  W]  (*  *V)W}  -  6((k*V)W,W(k*V)  Wj  -  6  ( [  ( »c  •  V )  W]  2,  W  ) 

-  18  ( (ic«V)  W  ,  W(k*  V)  w)  -  18  (W(tc*  V)  W,  (ic  *V)W)  . 

**  X 

Therefore, 

3t<W3  +  |  (k-VW)2>  -  3X<|  W4  +  12WtK*VW)2 
+  |  I(k-V)?W]2>  +  B<f [3W2-(k*V)2W]> 

*  ((k*V)  W,0cfc-7)W)  -  3((>c*V)2W,  (r*V)  (<X-V)W) 

-  ((k*V)W,  6VJ  (k«V)W  +  eW(kv»V)w)  (VII. 9) 

A  X 

=  ((K*V)W,  (<t«V)W  +  3(k«V)2(kx*V)w)  • 

-  6W  (k» V)  W  -  6W(k  »V)  W  . 

A  A 

Thus,  the  third  null  space  condition  will  yield  the 
third  averaged  conservation  law  provided  the  expression 

((K.«V)  W,  Oc.‘V)W+  3(k*V)  2(k  *V)  W-  6W  (k  •  V)  W-  6W(k  *V)  W)  (VII.  10) 

t  A  A  A 


vanishes.  To  reach  this  point,  we  merely  integrated  the  null 
space  condition  by  parts  several  times.  Now  we  employ  the 
first  N  modulation  equations 


-S 


to  replace  expression  (VII. 10)  by 


( (*c- V)  W,  (w  «V)  W  +  3(k*7)2(k*7)W  -  6W  (**7)W-  6W(*  *7)w) 

X  XX  X 


(VII. 10') 


To  show  this  expression  vanishes,  we  use  the  equation  satisfied 
by  the  N  phase  wave  s 


(w-V)W  -  6W(k«V)W  +  (<*7)JW  =  0  . 


Differentiating  this  equation  with  respect  to  x  yields 


(<*>  •  v )  W  -  6WU*7)WV  -  6W  (  V)  W  +  (k*V)jw 

■X  X 


-  (0)  .V)W  -  6W(k  *V)W  +  3(K‘7)4(ic  «7)W  =  0  . 

XX  X 

Finally,  take  the  inner  product  with  (k*7)W: 

(  (  K*  V )  W,  (a) •  7 )  W  -  6W (  <•  7 )  W  -  6W  (**7)W  +  (  K*7)3W 

X  XX  X 

+  (<*>  *7)W-  6W(  K  *7)W  +  3  (  k*7)  2  (k  *7)w)  -  0 
XX  X 


That  is. 


((k*V)Wx,  (o)»7)W-  6W (  k*  V) W  +  (k>7)?W) 

+  ((k'V)W,  (w  -V)W+  3  (  k*7)  2  (k  »7)W  -  6W  U‘7)W-  6W(*  *7)w) 

X  XX  X 


The  first  term  vanishes  by  the  fact  that  W  solves  the  N  phase 
equation;  thus,  we  arrive  at  the  vanishing  of  (VII. 10') I 


*V 


In  summary,  we  have  shown  that  (i)  the  first  two  null 


space  equations  imply  the  first  two  averaged  conservation  laws, 
and  that  (ii)  the  third  null  space  equation,  together  with 
the  consistency  conditions  =  wx  •  imply  the  third 
averaged  conservation  law. 

Remark  (i) ♦  Presumably,  the  j*"*1  null  space  equation, 
together  with  the  consistency  conditions,  implies  the  j 
averaged  conservation  law.  To  prove  this  statement,  one 
needs  a  more  abstract  argument  than  the  explicit  calculation 
described  above.  As  yet,  we  have  not  succeeded. 

Remark  (ii) .  We  have  completed  a  first  step  toward 
this  argument.  Consider  (6Hj/6w# Vfc) .  Using  the  identity  [31] 


H.  = 


_  1 


2  j  +  1 


OH.,. 


we  obtain 


8  <  1 
Tt  2 j+1 


a”i+i  -  a  „ 
Jw  it  ) 


3H.  3H . 

=  W3-'  Wt)  +  JT1  Kt  ; 


where  =  H ^ ( W; k) .  Thus, 


oh  .  .  a 

'  wt^  “  2j+r  ?t 


0  6H.^. ,  3H . 

<f  3.+1  >  -  _  —1  £ 

5t  5q  |  3  k  t 

W 


(VII. 11) 


Equivalently,  it  may  be  better  to  use  a  formula  from  [  ] , 


-  V  4-  — 

7W  "  *X  dx 


V 

r  ^X  «T*  _  VX  6¥~1 

(J  If  W"  +  2?  SWJ  ' 


(VII. 12) 


where  Yx  =  3^  ,  -  3  ▼  . 


. ..  .  ■.-.v.W •  .v.v. 


VIII.  An  Invariant  Form  of  the  Modulation  Equations 

In  this  section,  we  assume  the  correct  modulation 
equations  are  of  the  form 


3  k  .  =  3  a) .  j  =  1,2,...,N 

t  3  x  3  ' 


(VIII. la, b) 


3t  <H'j>  +  3x  (  Xj>  +  P<Gj>  j  =  1/2»-*wN+1 


We  summarize  some  results  of  [3]  in  order  to  emphasize  that 
these  modulation  equations  are  indeed  tractible. 


In  [3] ,  we  establish  the  following 


Theorem 


The  modulation  equations  (VIII. 1)  admit  an 


equivalent  representation  in  terms  of  the  differentials 
ni'fi2  : 


3^-12  3  «2  +  8  dF  -  0  , 


(VIII. 2) 


where  F  is  a  meromorphic  function  on  the  Riemann  surface 
R  of  the  form 


F  (X)  = 


A  ( 1 ) 

RW 


A  (X) 


IX 

*  s  “j  xj 
1  =  0 


(VIII. 3a) 


The  coefficients  a.  are  fixed  in  terms  of  <G.>  by  the 

3  3 

linear  system 


3  r  k  (-2N 


5+T  * Gj+i>'i-°'1'2' 


(VIII. 3b) 


Here 


-  £ 


r2k 

pk  5  • 


'  n  (1-Apr) 
1=0  * 


Representation  (VIII. 2)  is  fundamental.  It  contains 
alternative  representations  and  quickly  shows  they  are  equivalent. 
For  example,  the  most  useful  mathematical  form  is  an  immediate 


Corollary  (Riemann  Invariant  Form  of  the  Modulation  Equations) . 


By  evaluating  the  invariant  representation  -  9xfi2  +  &dF=0 

at  the  branch  points,  one  obtains 


(jn 

)tAA  +  S  (A)  3xA£  =  3  n+T - - -  »  £=0, 1,...,2N,  (VIII. 4a) 

L  J  S  D:x^_1 


where  the  characteristic  speed  s  ^  (A)  is  given  by 


38 


sU)  (X)  = 


-12  U  e3  x* 


(VIII. 4b) 


W Tj 

§vr 


A  fully  nonlinear  modulation  theory  cannot  be  simpler  than 

the  Riemann  invariant  form  (VIII. 4a).  We  emphasize  the  "internal 

perturbations"  described  in  the  introduction  cause  the  parameters 

ID 

to  modulate  with  the  characteristic  speeds  s  (X).  The 
external  perturbation  f  provides  the  right  hand  side  of  the 
modulation  equations. 


.  »  »  ■  4  *  *  *  V 

*  «  *  .  *  *  *  *  *  *  *  .  *  •  ▼  -  '  »  •  ’.«  4  ‘  V*.  *  '  ,  , 


Recently/  in  some  very  interesting  mathematical  work 
[19/20,10]  weak  limits  of  solutions  of  nonlinear  pde's 
have  been  described  by  a  measure.  When  the  nonlinear  pde 
is  dissipation  dominated,  as  in  Burgers'  equation, 
the  measure  is  simply  a  Dirac  measure  supported  on  the  shock. 
When  oscillations  persist  as  in  the  small  dispersion  limit 
of  KdV,  the  weak  limit  is  more  interesting.  As  yet,  the 
measure  has  not  been  rigorously  characterized. 

Our  purpose  in  this  final  section  is  to  calculate, 
with  heuristic  reasoning,  the  measure  which  describes  the 
weak  limit  as  e  -*•  0  of  problem  (II. la),  in  a  region  of 
space-time  where  the  solution  is  described  by  a  modulating 
N-phase  wave.  Since  modulation  theory  constructs  the 
solution  for  small  but  finite  e,  our  representation  certainly 
contains  enough  information  to  calculate  the  measure  very 
explicitly. 

Consider  the  solution  of 

Ut  -  6UUX  +  e2Uxxx  +  B  f(U,eUx,  ...)  =  0  ,  (I 

as  constructed  in  this  paper.  Namely,  in  a  region  S  of  space 
£ 

time,  let  U  be  the  modulating  N  phase  wave 


9  •  •  •  9 


e  re  (x,t) 

U  (x,t)  =  WNj-± - 


0N(x,t)  *  1 

- ^ -  ;  t(x,t)  j 

+  e  U(1)  ;  x,t  )+..., 


with  the  modulation  of  the  parameters  t(x, t)  described  by 


n  =  dF 


(IX. 2) 


(IX. 3) 


Fix  some  t  and  a  spatial  interval  I  such  that  (t,I)  G  S  ; 
further,  let  4> :  ]R  -*■  3R  denote  any  test  function  with 
support  within  I.  Fix  f:  3R  -»•  3R  and  consider 


co 

[♦,  f  (U(E)  ( .  ,t) )  j  =  }  4>  (x,t)  f  (U(E)  (x,t) )  dx  , 


which  one  must  consider  in  order  to  describe  the  weak  limit 
f (Ue(.,t))  as  £  +  0.  We  compute: 


w 

(*,f(u(e> <.,tn)  -  |  *(x,t)  f(u(e) (x,t>)  ax 


-  I  { 

Xs-oo  J 


xa+A/2 


4>  (x,  t)  f  ( IT  '  (x ,  t) )  dx 


L“  x±-  A/2 
X.+A/2 

CO  f  1 


*  l  f  4>(x,t)  f  t(trx))l  dx 

~  i—  .  WO 


x^- A/2 


Xi+V2 

+  l  <f»(x.,t)  j  ^(t.x^jidx 

xi“A/2  (A  tiny) 


f-  ‘•'V'-T? 


ihl 


00  f  n  x  r  rK(x.  ,t)x-3(x.,t)  ^  . 

I  „  4  |r  J  ftwNl — 1 — s - j! 


*i-v* 


OD  A/2£  1 

X  ♦  (x .  ,t)  A  -  ^  |  f[wN(*^y  ^(t,x1»]  dyl 

i=_“  -A/2£  J 


*lVl" 


.,®N ;  ^(xi,t)]dNe| 


(e  +  0  and  ergodicity) 


w 

^  |  0  (x ,  t )  <f(WN;  £(x,t))>  dx  (A  -V  0) 


Thus,  we  compute 


lim  (❖, f (UE(.,t))) 
e+0 


=  (<t>,  <£>  J  V  <j»  ; 


that  is. 


f  (U  (x ,  t ) )  -T^g- 


<f(WN(-;  t(x,t))> 


(IX. 4 


f  flWM(0?  t(x,t))I  dN0. 
iit)N  iN  N 


The  measure  itself  can  be  characterized  very  explicitly  by 


using  the  y-coordinates  for  the  torus  (see  VII. 18  ): 


f (Ue(x,t)) 


y  j  •••  J  £[A(Xj  t) 


TT  (p.-v,) 

i>l  1  3 

N 

I  I  |R()i|;t(xft))  | 
i=l  1 


dy,  A  ...  dp. 


(By  translation,  we  can  remove  the  function  A(x,t)  from  the 
argument  of  f.) 

Formula  (IX. 5)  is  the  main  result  of  this  section. 

It  shows  that  the  weak  limit  as  e  -*•  0  can  be  characterized 
by  a  measure,  and  gives  an  explicit  formula  for  the  measure 


Notice  that  the  (x,t)  dependence  of  the  measure  is 

through  t(x,t)  which  satisfies  ft  =  dF.  This  is  precisely  the 
measure  used  in  [1]. 


Remark .  Our  calculation  is  limited  to  the  region  S 
where  the  solution  is  an  N  phase  wave.  It  is  not  uniformly 
applicable.  However,  we  believe  it  indicates  that  the 
general  weak  limit,  which  would  be  valid  for  all  space, 
should  be  characterized  by  the  quadratic  variational  problem 
of  [8,93. 


*  » 
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